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> ■ ABSTRACT 

\ Using evolutionary population synthesis we present high resolution (0.3 A) integrated 

5^ ■ spectral energy distributions from 3000 to 7000 A and absorption-line indices defined 

by the Lick Observatory image dissector scanner (referred to as the Lick/IDS) sys- 
tem, for an extensive set of instantaneous burst binary stellar populations with binary 
interactions. The ages of the populations are in the range 1 — 15 Gyr and the metal- 
\ licities are in the range 0.004 — 0.03. This high resolution synthesis results can satisfy 

■"■^ — ■ the needs of modern spectroscopic galaxy surveys, and are available on request. 

' By comparing the synthetic continuum of populations at high and low resolution 

we show that there is a good agreement for solar metallicity and a tolerable dis- 
Q ■ agreement for non-solar metallicity. The strength of the Balmer lines at high spectral 

| H \ resolution is greater than that at low resolution for all metallicities. The comparison 

C/3 . of Lick/IDS absorption-line indices at low and high resolution, both of which arc ob- 

^ ' tained by the fitting functions, shows that the discrepancies in all indices except for 

TiOi and Ti02 are insignificant for populations with Z — 0.004 and Z — 0.02. The 
high resolution Ca4227, Fe5015 and Mgb indices are redder than the corresponding 
/\ ' low resolution one for populations with Z = 0.01 and Z = 0.03, this effect lowers the 

derived age and metallicity of the population. The high resolution Mgi, Fe5709 and 
Fe5782 indices are bluer than those at low resolution, it raises the age and metallic- 
ity. The discrepancy in these six indices is greater for populations with Z = 0.03 in 
comparison to Z = 0.01. 

At high resolution we compare the Lick/IDS spectral absorption indices obtained 
by using the fitting functions with those measured directly from the synthetic spectra, 
and see that Ca4455, Fe4668, Mgb and Na D indices obtained by the use of the fitting 
functions are redder for all metallicities, Fe5709 is redder at Z = 0.03 and becomes to 
be bluer at Z — 0.01 and 0.004, and other indices are bluer for all metallicities than 
the corresponding values measured directly from the synthetic spectra. 
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1 INTRODUCTION 

In the previous papers iZhang et al ] |2004l.l2005l hereinafter 
Paper I, II), we took into account various known classes 
of binary stars in evolutionary population synthesis (EPS) 
models, presented integrated colours, integrated spectral en- 
ergy distributions (ISEDs) and 21 absorption-line indices 
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[following the definitions of lWorthev et al1ll994ll defined by 
the Lick Observatory image dissector scanner (referred to as 
the Lick/IDS) system for an extensive set of instantaneous- 
burst binary stellar populations (BSPs) with binary inter- 
actions, investigated the influences of binary interactions 
and model input parameters on the results, and found that 
the inclusion of binary interactions makes the integrated 
U - B, B - V, V - R and R - I colours and the 21 Lick/IDS 
spectral indices substantially bluer. In Papers I and II the 
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corrected BaSeL-2 .0 (i.e., non-calibrated BaSeL-2.2) stellar 
spectra library of iLeieune. Cuisinier fc Buserl (Il997l Il998l) 
was used and low spectral resolution ISEDs (10 A in the 
ultraviolet and 20 A in the visible) were given. 

In recent years, some spectroscopic galaxy surveys at 
intermediate spectral resolution (such as the Sloan Digital 
Sky Survey, SDSS) have been undertaken, and the data has 
been analyzed with EPS models to understand galaxy for- 
mation and evolution. Therefore, a EPS model that can pre- 
dict the ISEDs at intermediate and high spectral resolution 
is required. The high spectral resolution of ISEDs can be 
used to predict the strength of numerous weak absorption 
lines and the evolution of t he profiles of the strongest line s 
over a wide range of ages jGonzalez Delgado et al]l2005l) . 
And, the higher t he spectral reso lution, the higher the con- 
straining power dVazdekislligggll . Therefore, in this paper 
we attempt to include high resolution spectral information 
in our EPS models. 

So far, the majority of EPS studies used low resolution 
st ellar spectral lib rary as an ingredient. In t he EP S models 
of lVazdekisI Jl999l) and iBruzual fc Charlol] J2003h interme- 



diate spectral resolution ISEDs were provided: the former 
used a subsample of ~ 500 stars from the original Ijoned 
1^93) empirical stellar library and presented the ISEDs of 
single stellar populations (SSPs) with a spectral resolution 
of ~ 1.8 A in two reduced spectral regions: 3856 — 4476 
and 4795 - 5465 A, the latter used the STELIB library 
jLe Borgne et aljEo OSl and presented 3 A resolution models 
covering 3200 - 9500 A with their code GALAXEV. In ad- 
dition, EPS models also ca n use the intermediate-resolutio n 
empirical stellar library of ISanchez-Bl azauez et all i2003l) . 
which comprises ~ 1100 stars and covers the spectral range 
3500 - 7500 A with a resolution of ~ 2 A. 

With the emergence of new-gcncration empirical 
the ELODIE data base of Prugnicl & Soubiran 
and 
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2004f) stellar spectral libraries, EPS models can present 
the high-resolution IS EDs and spectral absorption indices. 
iBertone et all ^l2003b^ used their library and explored the 
spectral properties of SSPs in the optical (4880— 5390 A) 
and uhraviolet regions (2275-2575 A and 2690-2960 A) 
at R = 500 000 and 50 00, respective l y. iTantalo et ail 
used the library of iMunari et al.l ^2005^ and pre- 
sented the Lick absorption-line in dices of SSPs from their 
1- A re solution spectra. Recently, iGonzalez Delgado et alJ 
i2005fl used their library and presented high resolution 
theoretical spectra (~ 0.3 A) of SSPs in the 3000 - 7000 A 
range. 

In the few EPS studies to date that have presented in- 
termediate and high spectral resolution results, binary in- 
teractions are neglected. Binary stars paly an important 
role in determining the overall appea rance of any realistic 
stellar population ( IXin fc DendbOOSl . and Papers I and II 
for an extensive discussion), therefore, in this paper we in- 
tend to use the high resolution stellar spectral library in our 
EPS study of BSPs and present their ISEDs and Lick/IDS 
absorpt ion-feature indices. Th is set of indices includes 21 in- 
dices of IWorthev et all il994ll and 4 Balmer indices defined 



bv lWorthev fc Ottavianil (^g^)- In an instantaneous burst 
BSP all stars are assumed to be born in binaries and born 
at the same time, and 1 x l(f binary systems are comprised. 

The outline of the paper is as follows: we describe our 
EPS models and algorithm in Section 2; our results are pre- 
sented in Section 3, and then finally, in Section 4, we give 
our conclusions. 



2 MODEL DESCRIPTION 

2.1 Initialization of the binary population 

We first need to construct the instantaneous burst BSPs. 
A Monte Carlo process is used to produce a population of 
1 X 10® binary systems, and the initial state of each binary 
satisfies the following input distributions: 

(i) the initial mass function (IMF) of the primaries, which 
gives the relative number of the primaries in the mass range 
mi — > mi -I- dmi. In this study the initial mass of the 
primary is chosen fr om the ap proximation to the IMF of 
iMiller fc Sc'alol ^1979^ as given bv lEggleton. Fitchett fc Tout! 
g,98j, hereinafter EFT), 



mi — 



0.19X 



(1 - X)0-75 + 0.032(1 - X)0-25 



(1) 



where X is a random variable uniformly distributed in the 
range [0,1], and mi is the primary mass in units of M©. 

(ii) the initial secondary-mass distribution, which is as- 
sumed to be correlated with the initial primary-mass distri- 
bution in this study. So, it depends on the initial primary- 
mass (as set by equation and the initial mass ratio, q, 
distri bution, which is assumed to be a uniform fo rm (EFT 
1989; iMazeh et al]ll992l: iGoldberg fc Mazehlll994l) . 



n{q) 



< g < 1, 



(2) 



where q — m2 /mi , m2 is the secondary mass in units of Mq . 

(iii) the distribution of orbital separations (or periods) . It 
is taken as constant in log a (where a is the separation) for 
wide binaries and falls off smoothly at close separations: 



an{a) = | 



where as, 



Uscpia/ao)''' , 



a < ao, 

ao < a < ai, 



(3) 



0.070, ao = 10 Re, ai = 5.75 x lO^R© and 



7/> « 1.2. 

(iv) the eccentricity distribution. A uniform form is as- 
sumed: 



where X is a random variable, as in Eq. 



(4) 



After the initial state (the masses of the component 
stars, mi and m2, the separation a and eccentricity e of the 
orbit) of a binary system in a BSP is set, we also need to set 
the lower and upper mass cut-offs mi and niu to the mass 
distributions and assign a metallicity Z to the stars: mi and 
rtiu are set as 0.1 and 100 Mq, respectively, Z = 0.004, 0.01, 
0.02, 0.03. The relative age r of the BSP is assigned within 
the range of 1 — 15 Gyr. 
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Figure 1. Number of binary systems experiencing RLOF as a 
function of time for solar-metallicity BSPs of 1 X 10® binaries. 
Solid, dashed, dot-dashed, dotted, thick-dashed and thick-solid 
lines represent the 1st, 2nd, 3rd, 4th, 5th and all RLOF, respec- 
tively. 



2.2 Input physics, parameters and algorithm 

W e use the rapid binary sta r evolution (BSE) algorithm 
of iHurlev. Tout fc Polsl (|2002) to evolve each binary in the 
BSP to an age of r, which gives us evolutionary param- 
eters such as the stellar luminosity L, effective tempera- 
ture Teff, radius R, current mass m and the ratio of ra- 
dius to Roche-lobe radius R/Rl for the component stars. 
Next we use the high resolution HRES stellar spectral li- 
brary of iGonzalez Deleado et alJ ll2005l) to transform the 
evolutionary parameters to stellar flux by interpolating the 
flux grid in the logTeff— logg— [Fe/H] plane, and then use 
equation (jSJ to obtain the monochromatic flux for an in- 
stantaneous BSP of a particular age and metallicity. For the 
spectral absorption-line indices in the Lick/IDS system, we 
obtain them by two methods: (i) measure them directly from 
the synthetic f l ux, (ii ) use the empirical fitting function s of 
IWorthev et"ai] |l993) and IWorthev fc Ottavianil lligQTt) to 
assign the Lick/IDS absorption indices to a star with a set 
of given evolutionary parameters, and then use equations 
®-l|Hll to der ive the indices for a BSP. 

The integrated monochromatic flux of a BSP is defined 

as 



E 

k=l 



(5) 



where fx is the SED of the fcth star. 

The integrated absorption feature index of the Lick/IDS 
system is a flux-weighted one. For the ith atomic absorption 
line, it is expressed in equivalent width {W, in angstroms). 



Wi,r,Z 



Y.l=i hex 



(6) 



where Wi is the equivalent width of the ith index of the 
fcth star, and fi,c\ is the continuum flux at the midpoint of 
the ith 'feature' passband. The local continuum for the i-th 
index is the run of flux defined by drawing a straight line 
from the midpoint of the blue pseudocontinuum level to the 
midpoint of the red pseudocontinuum level, the flux at the 
midpoint of pseudocontinuum is a average one and obtained 
by 



fi,P 



Ai 2 — Ai 1 



(7) 



where f\ is the stellar flux, as in Eq. 0, Ai,i and Ai,2 are 
the wavelength limits of the i-th pseudocontinuum sideband. 
For the ith molecular line, the feature index is expressed in 
magnitude. 



a. 



= -2.5 lot 



10 



-0.4ci 



fi,CX 



CA 



(8) 



where a is the magnitude of the ith index of the fcth star. 

Detailed descriptions of the BSE package of 
iHurig^^^^ I^Q P^^ ^^" ^ the empirical fitting functions 
of lworthgy^tall il99J) have been presented in Papers 
I and II, here we only mention several important input 
parameters required in the BSE code: the efficiency of 
common envelope ejection acE is taken as 1.0, the Reimers 
wind mass-loss coefficient is set constant at 0.3, and the 
tidal enhancement parameter B — 0.0. 

Following, we present a brief description of t he HR ES 
stellar spectral library of Gonzalez Dclgado ct ah] ll2005l) . A 
full discussion is given in lMartins et alJ l|2005j). This library 
includes the synthetic stellar spectra from 3000 to 7000 A 
with a final spectral sampling of 0.3 A. The spectra span a 
range of effective temperature from 3000 to 55000K, with 
variable steps from 500 to 2500K, and a surface gravity 
logg = —0.5 to 5.5 with dex steps of 0.25 and 0.5. For 
each temperature, the minimum gravity is set by the Ed- 
dington limit. The library covers several metallicities: twice 
solar, solar, half solar and 1/10 solar. Solar abundance ra- 
tios for all the elements, and a helium abundance of He/H 
= 0.1 by number are assumed. This high spectral reso- 
lution stellar library is based on Kuru cz local therm ody- 
namic equilibrium (LTE l atmospheres (lKuruczlll^93^ and 
the program SYNSPEC fHuben v. Lanz fc Jeffervl |l995'^ for 
stars wi th 8000 < T^g < 2700 K, the program SPEC- 
TRUM jGrav fc Corballvl 11994) and Kurucz atmospheres 
for stars with 4750 < T^a < 7750 K, non-LTE line-blanketed 
models for hot (27500 < T^a < 55000 K, Lanz fc Hubc^ 
2003), and PHOENIX LTE Une -blanketed models for cool 
stars (3000 < T^ s < 4500 K, iHauschildt fc BaronI 119991 : 
lAllard et al .1200 j) . This library is composed of 1650 spectra. 

By adopting the above set of input distributions and 
parameters, in Fig. we present the number of binary sys- 
tems experiencing Roche lobe overfiow (RLOF) as a func- 
tion of time for solar-metallicity BSPs of 1 x 10® binaries. 
Because some binary systems would experience RLOF for 
several times, in Fig. Q solid, dashed, dot-dashed, dotted 
and thick-dashed and thick-solid lines represent the 1st, 2nd, 
3rd, 4th, 5th and all RLOF, respectively. From Fig.^we see 
that RLOF mainly happens at early age, and the number 
decreases with time. During the past 15 Gyr ~ 11.6 per cent 
of the binaries would experience RLOF. 
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Figure 2. The high spectral resolution integrated spectral energy 
distributions as a function of age for solar-metallicity BSPs. From 
top to bottom, the ages r = 1, 2, 4, 7, 13 Gyr, respectively. 
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Figure 4. Comparison of the high resolution (solid line) with 
low resolution (dots) integrated spectral energy distributions for 
solar-metallicity BSPs at ages r = 1, 2, 4, 7 and 13 Gyr (from top 
to bottom, respectively). 

3 RESULTS 

3.1 The integrated spectral energy distribution 

In this part we present ISEDs covering the range 3000 — 
7000 A with a resolution of 0.3 A for instantaneous burst 
BSPs with binary interactions over a large range of age and 
metallicity: 1 < r < 15 Gyr and -1.3 < [Fe/H] < +0.2. For 
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Figure 5. Comparison between the high (solid line) and low res- 
olution (dots) integrated spectral energy distributions for r = 
1 Gyr BSPs with Z = 0.004, 0.01, 0.02 and 0.03 (from top to bot- 
tom, respectively). For the sake of clarity the fluxes at metallic- 
ities of 0.01, 0.02 and 0.03 are shifted downwards by an amount 
of 1.5. 



each model, a total of 1.0 x 10 binaries are evolved according 
to the algorithm given in the previous section. The full set 
ISEDs are available on request from the authors. 

In Fig. 121 we give the high spectral resolution ISED evo- 
lution for solar-metallicity BSPs at ages t — 1, 2, 4, 7 and 13 
Gyr in the wavelength range of 3000 - 7000 A. It shows that 
the continuum tends to be redder with increasing the age of 
the BSP, the variation in the shape of the continuum with 
time is significant at early age, and is almost constant at 
intermediate and late ages. Additionally, Fig. |21 shows that 
the strength of the Balmer lines decreases with the age of 
the BSP. 

The metallicity effect on ISED evolution of BSPs is 
shown in the top panel of Fig.|3 which presents the ISEDs of 
young (r = 1 Gyr) and intermediate-age {t — 7 Gyr) BSPs 
at four different metallicities: Z = 0.004, 0.01, 0.02, 0.03. The 
metallicity effect is also exhibited in the slope of the contin- 
uum and in the strength of the metallic lines. The continuum 
tends to be redder with increasing metallicity and metallic 
lines to be stronger. The bottom panel of Fig. |H] gives the 
evolution of Bamler lines with metallicity, it shows that their 
strength tends to be weaker when increasing metallicity. 

In Fig.^Jwe compare the ISEDs at high and low spectral 
resolution (10 A in the ultraviolet and 20 A in the visible) 
for solar-metallicity BSPs at ages r = 1, 2, 4, 7 and 13 Gyr. 
These low resolution ISEDs are from Model A of Paper II, 
which includes binary interactions and adopts the same in- 
put parameters and distributions as this paper except for us- 
ing the corrected BaSeL-2. (i.e., non-cal i brated BaSeL-2.2) 
stellar spectral library of iLeieune et alJ il997l . llQQSl) . Fig. 
21 shows that there is a good agreement in the continuum 
for solar-metallicity BSPs between two studies. Comparison 




of the continuum for BSPs with non-solar metaUicity shows 
that there is a difference existed at longer wavelengths for 
BSPs at ages r = 1 and 2 Gyr. In Fig^we give the ISEDs at 
high and low resolution for r = 1 Gyr ESP with four metal- 
licities. It shows that at longer wavelengths the high resolu- 
tion continuum is lower (redder) for BSPs with Z — 0.004, 
while greater (bluer) for BSPs with Z = 0.01 and 0.03 than 
that at low resolution. 



3.2 Lick spectral absorption feature indices 

At high spectral resolution we use two methods to ob- 
tain 25 Lick/IDS sp ectr al absorption indices define d by 
IWorthev et all Jl994) and lWorthev fc Qttavianl JlOQl) : (1) 
directly compute them from the high-resolution synthetic 
spectra, (2) obtain them by usi ng the empirical fitting func- 
tions of IWorthev et alJ and IWorthev fc Ottaviaiiil 
Jl997fl . 

All the following results are for BSPs of lO'' binaries. 
To verify that our results are stable for 10^ binaries, we 



perform several simulations, the number of binaries in each 
simulation is 10*, 4 x 10*, 1.6 x 10^ 3.6 x 10^ and 6.4 x 10^ 
respectively. As an example, in Fig.|S|we plot the evolution 
of Mgi index as a function of the number of binaries in the 
simulation, we see that the Mgi index fluctuates around the 
value for 10® binaries when the number of binaries is not 
so enough, with the increase of the number of binaries the 
fluctuation would decrease. For 10* binaries the discrepancy 
reaches to ~ 0.05 A at r = 12 Gyr, for 3.6 x 10^ binaries Mgi 
agrees with the value for 10® binaries within a small range, 
and for 6.4 x 10^ binaries the discrepancy is insignificant. 
Therefore the results for 10® binaries in the simulation are 
stable. 



3.2.1 Results by the use of the fitting functions 

U sing the high re s olutio n sp ectra and the fitting function s 
of IWorthev et alJ ^1994^ and IWorthev fc Ottaviaiiil Jl997D . 
we obtain 25 Lick/IDS spectral absorption indices for the 
BSPs with binary interactions with four metallicities from 
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Table 1. Lick/IDS spectral absorption feature indices of BSPs, derived by using the empirical fitting functions (Section 3.2.1). 
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Figure 7. Evolution of absorption indices in the Lick/IDS system obtained by using the fitting fun ctions for BSPs of various metal licity. 
The symbols linked by a line denote the indices obtained b y using the hi g h resolution spectra of iGonzalez Delgado et al.l j200^, and 
those without a line are from the low resolution spectra of lLeieune et all ll997l [l998l) . Different symbols are for different metallicity, 
from top to bottom, the metallicity Z is 0.03, 0.02, 0.01 and 0.004, respectively. 



Z = 0.004 to Z = 0.03. The results are presented in Table 
Q We compare them with those obtained by Model A of 
Paper II, which uses low resolution stellar spectral library 
of and the fitting functions, and find that the evolution of 
TiOi and Ti02 indices is not smooth for both two studies, 
the discrepancy in the rest 23 Lick/IDS spectral indices is 
small for metallicity Z = 0.004 and Z — 0.02, the synthetic 
Ca4227 (index 3), Fe5015 (10), Mgi (11), Mgb (13), Fe5709 

(17) and Fe5782 (18) indices show much stronger discrep- 
ancy for metallicity Z = 0.01 and Z — 0.03 . 

In Fig. Qwe show the comparison of Ca4227 (index 3), 
Fe5015 (10), Mgi (11), Mgb (13), Fe5709 (17) and Fe5782 

(18) indices between two studies. It shows that the high 
spectral resolution Ca4227 (3), Fe5015 (10) and Mgb (13) 
indices are greater, Mgi (11), Fe5709 (17) and Fe5782 (18) 
indices are less than those at low resolution for BSPs with 
Z = 0.03 and Z = 0.01, the discrepancy is greater for BSPs 
with Z = 0.03 in comparison to Z = 0.01. This discrepancy 
in the Lick/IDS spectral indices will directly influence the 
derived age and metallicity of a particular population, it 
will lower the age and metallicity if using high resolution 
Ca4227 (3), Fe5015 (10) and Mgb (13) indices, and raise the 
age and metallicity if using high resolution Mgi (11), Fe5709 
(17) and Fe5782 (18) indices. 

To discuss whether the difference of the coverage in the 
log(?— logTeff plane between HRES and BaSeL-2.0 libraries 
would cause the discrepancies in the Lick/IDS spectral ab- 
sorption indices, in Figs. |H| and |5| we plot the theoretical 



isochrones of the bluest (i.e., the youngest [r — 1 Gyr] and 
lowest metallicity [Z = 0.004]) and the reddest (the oldest 
[r = 15 Gyr] and highest metallicity Z = 0.03) BSPs. In or- 
der to transform the evolutionary parameters of stars along 
the isochrones to stellar flux for Z = 0.004 BSPs in Fig. 
|H| we need to make interpolation in the flux grids between 
[Fe/H]= —1.0 and —0.5 because HRES spectral library only 
contains [Fe/H] = -1.0,-0.5,0.0 and 0.3. While, the grid 
coverage in the logg— logTcff plane between [Fe/H] = — 1.0 
and —0.5 is appreciably different, and the difference mainly 
concentrates in the low-temperature range. Within this low 
temperature range the denser stars (such as, low-mass MS 
stars) only give small contribution to the ISEDs (see Figs. 
3 and 4 of Paper H), so we only interest in the difference of 
two libraries in the low-density range, in this range the grid 
coverage at [Fe/H] = —0.5 is smaller than that at [Fe/H] 
= —1.0, so in Fig. |51we plot the boundary of HRES flux 
grid coverage at [Fe/H] — —0.5. Also we plot the bound- 
ary of BaSeL-2.0 flux grid coverage at [Fe/H] — —0.5, for 
this library the grid coverage in the logp— logTcff plane for 
[Fe/H] = —1.0 and —0.5 is same in the low-temperature and 
low-density range. Similarly, for Z — 0.03 BSPs in Fig. 1^1 
we plot the boundary of HRES flux grid coverage at [Fe/H] 
= 0.3 and BaSeL-2.0 library at [Fe/H] = 0.2. 

Comparison of the coverage in the logg— logTeff plane 
for two libraries in Figs.|S]and|5]shows HRES library extends 
to the higher temperature, in the logTcff > 4.43 range covers 
less gravity than BaSeL-2.0 library, while the latter covers 
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Figure 6. Sensibility of Mgi index to tlie number of binary sys- 
tems in the BSPs. 



lower temperature range. Although in the high temperature 
range the grid coverage of two libraries exists discrepancies, 
both grid coverage are enough to the BSPs we considered. 
In the lower temperature and lower gravity range, HRES 
library does not cover all cool AGB stars, it would cause 
the discrepancies in the Lick/IDS absorption line indices. 
Additionally, the differences in the spectral shape caused by 
the resolution also would cause to the discrepancies of the 
above several Lick/IDS absorption line indices. 



3.2.2 Results computed from the high-resolution spectra 

To compare theoretical results with observations, the 
Lick/IDS absorption indices also are measured directly 
from the high-resolution ISED s usi ng the passband defi- 
nitions of IWorthev et alJ (^3) and IWorthev fc Ottavianil 
l|l997|) . and are presented in Table |5| In Fig. 1101 we com- 
pare them (except for TiOi and Ti02) with tho s e obt ained 
by using the fitting function s of IWorthev et afl I^O^) and 
IWorthev fc Ottavianil 1^93) (generated from Tabled, and 
find that Ca4455 (index 6, almost equal for Z — 0.02), 
Fe4668 (index 8, almost equal for Z = 0.004), Mgb (13) and 
Na D (19) are bluer than the corresponding ones obtained 
by using the fitting function for all metallicities. For Fe5709 
(17) the discrepancy becomes from negative {Z — 0.03) to 
positive {Z =0.01 and 0.004). Other indices are redder than 
the values from Table Q for all metallicities. 

Comparing the discrepancies in the Lick/IDS spectra 
indices caused by the difference in the metallicity, we find 
that the discrepancies in Ca4455 (index 6), Mg2 (12), Mgt 
(13) and Na D (19) indices introduced by using the difi'erent 
computation method are smaller. 
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Figure 8. Coverage of HRES (solid line) and BaSeL-2.0 (dashed 
line) stellar library grids for [Fe/H]= —0.5. For the sake of clar- 
ity BaSeL-2.0 library is shifted upwards by an amount of 0.1. 
Over-plotted points represent the theoretical isochrone for in- 
stantaneous burst BSPs with Z = 0.004 and r = 1 Gyr. Only 
200,000 binary systems are plotted and those MS stars with mass 
A/<0.7Mq are removed. 




Figure 9. Similar to Fig. but [Fe/H]= 0.3 for HRES and 
[Fe/H]= 0.2 for BaSeL-2.0 libraries, Z = 0.03 and r = 15 Gyr 
for the BSPs. 
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Table 2. Lick/IDS spectral absorption feature indices of BSPs, as measured directly from the high-resolution ISEDs (Section 3.2.2). 
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Figure 10. Absorption indices in the Lick/IDS system computed directly from the high-resolution ISEDs (symbols + line) and by using 
the fitting functions (only symbols) for BSPs of various mctallicity. The symbols have the same meaning as in Fig. |7| 
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4 SUMMARY AND CONCLUSIONS 

We have simulated realistic stellar populations composed of 
100 per cent binaries by producing 1 x 10^ binary systems 
using a Monte Carlo technique. Using the EPS method we 
present the high spectral resolution (0.3 A) ISEDs over the 
wavelength range 3000-7000 A and 25 Lick/IDS spectral ab- 
sorption indices, for an extensive set of instantaneous burst 
BSPs with binary interactions over a large range of age and 
metallicity: 1 < r < 15 Gyr and 0.004 < Z < 0.03. This set 
of high spectral resolution ISEDs fully satisfies the need of 
studying the formation and evolution of galaxy by analyzing 
the data of modern spectroscopic galaxy surveys, and can 
be available on request. 

By comparing the synthetic continuum at high and low 
resolution, we show that there is a good agreement for BSPs 
with Z — 0.02 and smaller discrepancy for BSPs with non- 
solar metallicity. And, the comparison of the Lick/IDS spec- 
tral indices at low and high resolution, both of which are 
obtained by using the fitting functions, shows that the high 
resolution Ca4227, Fe5015 and Mgb indices are redder, Mgi, 
Fe5709 and Fe5782 indices are bluer than those at low res- 
olution for BSPs with Z = 0.01 and 0.03. The high spectral 
Ca4227, Fe5015 and Mgb indices act to lower the age and 
metallicity, and the high resolution Mgi , Fe5709 and Fe5782 
indices will raise the age and metallicity. 

At high resolution we compare the Lick/IDS spectral 
absorption indices obtained by using the fitting functions 
with those measured directly from the synthetic spectra, 
and find that Ca4455, Fe4668, Mgb and Na D indices are 
redder for all metallicities, Fe5709 is redder at Z = 0.03 and 
becomes to be bluer at Z =0.01 and 0.004, other indices 
are bluer for all metallicities than the corresponding values 
computed from the high-resolution ISEDs. 
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